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Recent developments in theory and phenomenology relevant to deep inelastic lepton 
, scattering are reviewed, concentrating on the following topics: predicted behaviour 

of non-singlet and polarized structure functions at small x\ theoretical studies of 
saturation and unitarity effects at small x in quarkonium scattering; renormalons 
and higher twist contributions; next-to-leading-order calculations of jet cross sec- 
tions; forward jet production as a probe of small- a: dynamics. 

, 1 Introduction 

o 



The aim of this introductory review is to highUght some theoretical and phe- 
nomenological work done in the past year or so which may be relevant to deep 
inelastic scattering (DIS). The selection of topics is of course subjective; it 
is intended to cover a range of physics somewhat superficially, relying on the 
Oh, talks to be presented in the working groups to fill in the details. 

In the important field of small- a; physics, I have chosen to concentrate on 
rather non-standard topics: the small-x behaviour of non-singlet and polarized 
structure functions, and the small- a: dynamics of (quark)onium. This is be- 
, cause there has been interesting new work in these areas, whereas in the main 

' arena of unpolarized singlet structure we are experiencing something of a lull 

until more precise calculations become available. 

On the large topic of diffractive DIS I also have little to say, although 
the study of multi-pomeron exchange in onium scattering is relevant. The 
main developments in this area are experimental at present and substantial 
theoretical progress may be expected to come later. 

I have included a rather speculative section on renormalons, because this 
is a fashionable topic throughout QCD at present, and it might just give some 
indication of where progress could be made on the difficult problem of higher- 
twist contributions. 



"Talk at DIS96, International Workshop on Deep Inelastic Scattering and Related Phe- 
nomena, Rome, April 1996. 
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Finally I discuss jet production, for which important new calculations have 
been completed, although the situation on their comparison with older calcu- 
lations and with experiment is still unclear. 



2 Non-Singlet and Polarized Structure Functions at Small x 



There is a lot of interesting physics associated with the behaviour of non-singiet 
and polarized structure functions, for example that connected with sum rules. 
The Gottfried sum 



1 + / dx(u - d) 





(1) 



reveals that the u and d distributions are not the same, while the Gross- 
Llewellyn-Smith and Bjorken sum rules 



dx{gl - g^) 
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provide precise determinations of a^. The Ellis- Jaffe sum rule for g\ alone 
suggested the "proton spin crisis". All these sum rules require assumptions 
about how to extrapolate non-singlet quantities towards a; = 0, so we need to 
understand just what QCD predicts on this point. This is an area where there 
has been recent progress in theory and phenomenology. 

First let us recall what happens in the unpolarized, singlet case. The be- 
haviour is dominated by ladder graphs with longitudinally polarized gluons on 
each side. Fig. ^a), which give a basic x~'^ dependence to the gluon distribu- 
tion g(a;)_and hence to the structure functions. In DGLAP (leading-log Q^) 
evolution^ this is modified by a relatively slowly- varying asymptotic factor, 



1 



■ exp 



TT 



In- 



(DGLAP). 



(4) 



The presence of two large logarithms for each power of ois is characteristic. 
Here one keeps only those terms in which one of them is a log . These come 
from the region of phase space in which the transverse momenta fcr of the 
exchanged gluons are ordered along the ladder. 

At small X we need to keep instead the terms with the largest number of 
factors of logx. The fcr-ordering breaks down, and there is the possibility of 
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Figure 1: (a) Gluon ladder and (b) quark ladder graphs, contributing to singlet and non- 
singlet structure functions respectively. 



double-logarithmic (as In^ a;)" terms. However, BFKLq showed that double 
logs of X cancel in the fully inclusive, unpolarized, singlet structure functions. 
Although kT is not ordered, it diffuses without violent fluctuations. This leads 
to the asymptotic behaviour 

g{x) ^ - exp (xsCA—\n -] ^ (BFKL). (5) 

X \ TT X J 

The leading singlet BFKL eigenvalue xs = 4 In 2 implies that 75 ~ 0.5 for 
as ~ 0.2. 

The observed behaviour is consistent with next-to-leading-order (NLO) 
DGLAP evolution, and with the leading-order BFKL form (||) for a reduced 
value of 75. It is expected that non- leading corrections to the BFKL prediction 
will be importanJ.. They are not yet known, but progress is being made on 
calculating thcmS 

2.1 Non-Singlet Structure Functions at Small x 

For a non-singlet quantity /(a;), such as {F2 — F2) /x, we have quarks instead of 
gluons at the sides of the ladder, as in Fig. ||(b), and leading log-Q^ evolution 
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now gives 



/(a;)^i.exp( j2C7F^ln^lni) (DGLAP). (6) 



However, in this case the leading log-x terms are double-logarithmici'i Using 
a Sudakov representation of the exchanged momenta, 

fcM = xpf" + yqf + kl^ , (7) 

where q = q -\- xp (so that cp = 0), the double logarithms come from the 
region in which both momentum fractions are strongly ordered, but in opposite 
directions: 

x < x' < x" < • • • , y > y' > y" > • • • . (8) 
The dominant region of transverse momenta now becomes 

kl/x > k%/x' > k"l/x" > ■ ■ ■ . (9) 

Thus there can be huge fluctuations along the ladder, and the breakdown of 
kx ordering is much more drastic than in the singlet case. This leads to the 
Kirschner-Lipatovlj (KL) type of small-x behaviour: 



x-^Ns (KL) (10) 



where = ^J^Cpo^Jr; ~ 0.4 for ~ 0.2. Note that Regge behaviour at 
small X would imply that 

/(x) ~ x^"" with ap ~ i (R), (11) 

which is similar. 




2.2 Polarized Structure Functions at Small x 

The structure function combination Agi = g\ — 5" appearing in the Bjorken 
sum rule depends on the odd-signature part of the forward virtual Compton 
amplitude (i.e. odd under s-u crossing). Therefore the expected Regge be- 
haviour is that associated with exchange of the odd-signature, isovector ai 
trajectory, 

Agi(x) - withaai^iO (R). (12) 
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Figure 2: Non-ladder contributions that cancel for even-signature structure functions but 

not for odd signature. 

Using the KL approach^ however, Bartels, Ermolaev and Ryskini (BER) find 
that the corresponding power 7^ J is actuaUy shghtly larger than that for even 
signature: 

Agi(x) - > a;-^^^,' ^ x~°-^ (BER), (13) 

and so the behaviour as a; — *■ is more singular. 

The reason for the difference between even and odd signature is that there 
are non-ladder contributions, like those in Fig. |^, which cancel in the former 
case but not in the latter. The even-signature cancellation follows from the fact 
that the system exchanged in the t-channel is a colour singlet. The surprising 
thing is that the non-ladder terms enhance the odd-signature contribution 
rather than reducing it. 

An investigatimi of the phenomenological consequences for Agi has been 
presented recentlyllj Fig. ^ shows the results of resumming the large logarithms 
of X as discussed above, compared with those of the standard next-to-leading 
order DGLAP evolution. A difficulty with the resummation is that the leading 
log- a; expressions do not conserve fermion number, i.e. the integrals of the 
resummed splitting functions do not vanish. One can fix this in various ad hoc 
ways which do not affect the small- a; behaviour. The most extreme (curve A 
in Fig. ^ is to add a (5-Xunction contribution at x = 1. Curve B represents a 
less drastic prescriptionKi We see that the effect in either case is not great at 
foreseeable values of x and . 
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Figure 3: Small-x evolution oi — g" . A and B denote different prescriptions for fermion 

number conservation (see text). 



The singlet contribution to gi at small x has also been considered0 Here 
gluon ladders like Fig. |l|(a) can contribute, but only one of the gluons on the 
sides can be longitudinal, and so the basic power behaviour is = x", the 

same as for quark ladders, Fig. |l|(b). Again, the odd signature of gi means 
that non-ladder contributions like those in Fig. ^ do not cancel and in fact 
enhance the structure function. The asymptotic form is 



gi{x) ~ exp Cs 



TT X 



(BER). 



(14) 



The gluonic contribution alone would give Cs 
1.2. Mixing with quark exchange reduces the leading eigenvalue slightly: 



4 and hence 7^ 



37^ 



Cs 



3.5 



1.0 



(15) 



implying that 51(2;) ~ x^^. Further experimental data on the small-x be- 
haviour |Of gi are needed to test this surprising prediction. The theoretical 
analysis til suggests other remarkable features, such as a possible change in the 
sign of the polarized gluon distribution between large and small x. 
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Figure 4: Development of dipole chain in onium wave function. 

3 Small-a; Saturation and Unitarity in Onium Onium Scattering 

Amongst the important new phenomena expected in deep ineJastic scattering 
at small x are parton saturation and unitarity correctionstj As the parton 
density increases, interactions between partons become more significant and 
limit the rate of growth of the density (saturation). In addition, multiple 
scattering effects become necessary in order to prevent cross sections from 
exceeding unitarity limits. These phenomena are difficult to predict reliably in 
the case of ep scattering because the proton is a large object whose structure 
cannot be computed perturbatively. To study them within the domain of 
perturbation theory we need to consider a smaller system such as a heavy 
quark-antiquark bound state. For a sufhciently large quark mass mg, the size 
of such a state is i? ~ l/asTOg. Then the strong coupling at the relevant 
momentum scale ^ 1 /i?^ will be small enough for a perturbative expansion 
in as{q^) to be applicable. 

The theory of low-x phenomena in quarkonium (often simply called onium) 
has been worked out in aoma^detail over the last couple of years using the 
colour dipole formulation E^~E3 of QCD at small x. The dipole formulation is 
based on the fact that, to leading order in the number of colours, the light-cone 
wave function of an onium state can be represented by a chain of colour dipoles 
linking the heavy quark and antiquark (Fig. As the wave function is probed 
at smaller x, i.e. higher 'rapidity' y = ln(l/a:), more and more dipoles are likely 
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to be resolved. The evolution of the dipolc configurations with increasing y 
can be formulated as a branching process governed by the BFKL evolution 
equation. It turns out to be most convenient to describe the wave function 
in a (transverse position, rapidity) (r,y) representation. The mean density of 
dipoles of size c at a distance r from a primary onium {QQ) dipole of size h is 
found to be 

(16) 



n(&, c,r,y) - exp 
where 



ln2(16rV6c) 

(ap - l)y - ■ 



ky 



ap = 1 + 41n2CAas/7r (17) 
is the BFKL pomeron intercept and 

k^lAC{-i)CAash (18) 

is the corresponding diffusion constant. Thus as y increases the dipole density 
grows in the expected manner and dipoles diffuse slowly (logarithmically in r) 
from the onium source. The growth of dipole density leads to saturation and 
unitarity effects that can be investigated using the 'theoretical laboratory' of 
onium physics :[] 



3.1 Onium-Onium Scattering 

To study tie effects of unitarity it is simplest to consider onium-onium elastic 
scattering ll^^im at very high energy ^/s. The onia pass through each other 
with relative rapidity Y ^ ln(s/TOQ) at some relative impact parameter r. In 
the centre-of-mass frame, the dominant elastic process is two-gluon exchange 
between dipoles in the two onium wave functions evolved to rapidity y = Y/2. 
This corresponds to exchange of a single BFKL pomeron: the two evolved 
wave functions provide the upper and lower halves of a gluon ladder, which 
are joined in the middle by the two exchanged gluons. The two-gluon exchange 
interaction is short-ranged (it falls off like 1 /r'* at large transverse distances) , 
and so only those dipole that overlap contribute significantly to the scattering 
amplitude. 

Fig. m illustrates the situation for the two typical onium configurations 
shown at the top|| The impact parameter profile of the onc-pomeron ampli- 
tude, shown lower left, is obtained by moving one configuration over the other 

''Applica*i«nB-| of the dipole formulation to DIS phenomenology are also being 
invostigatedt3~a i—i 

'^This figure and tb« results that follow were obtained by Salamtj using his Monte Carlo 
program OEDIPUS Ej (Onium Evolution, Dipole Interaction, Perturbative Unitarization 
Software) . 
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Figure 5: Typical onium dipole configurations and their one-pomeron and unitarized inter- 
action amplitudes. 
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and computing the sum of the two-giuon exchange amphtudes between all 
pairs of dipoles at each relative position. Thus at impact parameters for which 
many dipoles overlap the one-pomeron amplitude violates the unitarity limit, 
as shown by the spikes exceeding unity. 

3.2 Unitarity and Saturation 

Schematically, the elastic ^-matrix element for single pomeron exchange has 
the form 

5=l-^P^Py^.y (19) 
7.7' 

where 7 and 7' represent dipole configurations of the two onia with probabilities 
Pj and P-yi respectively, and i/-y.^' is the one-pomeron amplitude for those 
configurations. To satisfy unitarity we need to take into account multiple 
pomeron exchange. For large dipole multiplicities the n-pomeron amplitude 
exponentiates, i.e. 

^ = i-EE^-^^r4:l (20) 

n 7,7' 



where 



and hence 



4:1 - J-Ji^ (21) 



5~E^'r^'r'exp(-/7,7') • (22) 



This ensures that the elastic amplitude never exceeds the unitarity limit, as 
illustrated for the particular configurations in Fig. |^ by the profile on the lower 
right. The resulting changes in the total and elastic cross sections as functions 
of rapidity are shown in Fig. ^. 

We see that the one-pomeron total and elastic cross sections have the 
expected energy dependences, s"^^^ and s^*^"^~^' respectively. The unitarized 
total cross section has a somewhat reduced effective power, while the elastic 
is more strongly affected. This is because elastic scattering is more dominated 
by smaller impact parameters, where large dipole densities are more probable. 

When one performs the sum over multiple pomeron exchanges first, and 
then the sum over dipole configurations, as in Eq. (^), it is obvious that a 
stable, unitary result is obtained. On the other hand one may evaluate the 



terms in Eq. (20) separately for each n, by summing over configurations fast. 
However, the multiple scattering series so obtained is strongly divergent£3^ll3 
The reason is that the configuration probability distributions and P^i have 
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Figure 6: Onium-onium total and elastic scattering cross sections. 



long (exponential) tails at high dipole multiplicities, which become dominant 
at high n but are strongly suppressed after summation over n. This suggests 
that for light hadrons also an expansion in terms of the number of pomerons 
exchanged may not be useful at very high energies. 

So far, we have discussed everything in the onium-onium centre-of-mass 
frame, where the two onia evolve through equal rapidity intervals, y = Y/2. We 
could equally well consider a 'fixed-target' type of frame, in which one onium 
evolves through y ^ Y and the other remains close to a pure QQ configuration. 
One can show that the above treatment is frame-independent at the level of 
one-pomeron (two-gluon) exchange, but the results for multiple exchange are 
frame dependent. The reason is that multiple scattering in the cm. frame 
looks like parton-parton interaction (saturation) in the wave function of the 
moving onium as seen from the rest frame of the other. Thus unitarity and 
saturation effects are intimately related, and one can learn about saturation by 
enforcing the condition that the simultaneous treatment of the two phenomena 
should be frame-independent EJ The reader is referred to the talk by Salamc3 
for further details. 
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4 Renormalons and Higher Twist 

An important thing to bear in mind when computing predictions of perturba- 
tive QCD is that the perturbation series is not expected to converge. There 
may already be indications of this in some quantities relevant to DIS: for ex- 
ample, the perturbative coefficients in the Gross-Llewellyn-Smith and Bjorken 
sum rules, Eqs. (J^ and (^), are growing rapidly. 

A known source of divergence of the perturbative expansion is the set of 
so-called renormalon graphsE^ illustrated in Fig. ^. Summing chains of vacuum 
polarization bubbles on the gluon line leads to a prediction of the form 

F = 5]c„a: (23) 

n 

where the coefficients Cn are factorially divergent at high orders: 

c„ ~ n\ {b/pY , (24) 

p being a rational number that depends on the particular quantity F being 
evaluated. Including n/ flavours of fermionic bubbles like those shown in Fig. ^ 
one finds h = —nf/Qn. To take_accQunt of gluonic bubbles as well, the 'naive 
non-Abelianization' procedure EJ~e3 is generally adopted: one evaluates the 
graphs with quark bubbles and then simply makes the replacement 

33 

nf^nj- — , (25) 
i.e. one uses for h the full 1-loop QCD /3-function coefficient h — (33 — 2nj)/127r. 



In fact the /3-function in QCD (unlike QED) cannot be associated with 
vacuum polarization graphs alone, and there is no set of graphs that cor- 
responds precisely to Eq. ( p3| ) with b = (33 — 2n/)/127r. Nevertheless the 
naive non-Abelianization procedure seems to work quite well in predicting the 
high-order coefficients in perturbation theory. An impressive example is the 
weak-jx)upling expansion of the average plaquette in quenched (n f = 0) lattice 
QCDEJ For this quantity the perturbative coefficients c„ in Eq. ( P3|) up to 
n = 8, suitably re-interpreted to allow for the diffepepce between the lattice 
and continuum definitions of a^, have been shown til to be dominated by a 
renormalon contribution of the form ( |2^ ) with b = 11/Att and p = 2. 

Although the series ( |2^ ) is divergent, it is not useless: it is expected to be 
an asymptotic expansion. By this we mean that the series truncated at order 
n = m is an approximation to the true value of F, with an error bounded by 
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Figure 7: A graph containing a "renormalon chain" 



the last term retained, CmCt™. Now the coefficients c„ satisfy 

-^^n- (26) 
c„-i p 

and thus the smallest term occurs at n = m ~ p/bag. Truncating the series at 
this point, we obtain the minimal ambiguity in the prediction of F, which is 

5F - m! m-" - ~ exp ( -7^ ) . (27) 



6a., 



Noting that 

1 



' 61n(Q2/A2) 
where is the relevant hard scattering scale, we see that 



(28) 



SF^[^\ . (29) 



Q 



We may distinguish two cases: 



1. p < 0. In QCD, this is an ultraviolet renormalon. The divergent series 
( p3| ) has alternating signs. It might appear from Eq. ( p9|) that this will 
lead to an ambiguity in F that grows with increasing Q^. In fact, how- 
ever, the series can be summed unambiguously in this case and there is 
no renormalon ambiguity. 

2. p > 0, which is called an infrared renormalon. All the terms of the diver- 
gent series have the same sign and it cannot be summed unambiguously. 
The ambiguity (E9) decreases like Q"^^ with increasing Q"^. 
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Figure 8: Power corrections to F2 (solid) and XF3 (dashed). Points are values deduced from 
BCDMS and SLAC data on F2 (see text). 



The above general discussion has some interesting possible applications to 
deep inelastic scattering. In the case of ultraviolet renormalons, various tech- 
niques can be used to sum theJivergent terms and thus improve the precision 
of the perturbative predictiono^Lj Infrared renormalons, on the other hand, 
imply that the perturbative prediction must be supplemented by additional 
non-perturbative information. In DIS, the power-suppressed nature of the 
ambiguity suggests that the extra information required concerns higher-twist 
contributions. In fact one can showtHI that similar renormalon graphs generate 
an ambiguity in the coefficient of the leading Q^^^ (twist 2p + 2) contribution, 
which cancels the power-suppressed ambiguity in the twist-2 prediction. 

Thus the presence of an infrared renormalon signals the need for a higher- 
twist contribution and tells us what its power dependence on must be. This 
suggests that we could go further and explore the possibility that the Bjorken-cc 
dependence of the renormalon also indicates the ^-dependence ofj-tlKiiigher- 
twist contribution. PhenomenoJogic.ally, this appears to be the caseE3^E3 Fig.|| 
shows, for example, predictionsE3'E3 based on this idea for the 1 / Q^-suppressed 
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corrections to the structure functions F2 and xF^, expressed in the form 



F{x,Q')~xq{x,Q') 1 + 



(30) 



where q is the appropriate combination of quariLdistributions. The form of the 
coefEcient D2 for F2 is similar to that deducfijdizl from BCDMsS and SLACB 
data, shown by the points. Recent results e3 on higher twist in xF^ also look 
similar in shape to the relevant curve (dashed) . It would clearly be of interest 
to include the type of higher-twist contributions suggested by this approach in 
global fits to deep inelastic scattering data. 



5 Jet Production 

The important topic of jet production in DIS has also seen recent theoretical 
and phenomenological progress. There are new results from complete next-to- 
leading-order (NLO) calculations of iet rates, which will be useful iar improving 
measurements of the strong couplingcll and the giuon distributioneS In addition 
there are new calculations on testing BFKL dynamics at small x with forward 
jet production as a trigger. 



5.1 Jet Cross Sections in Next-to- Leading Order 

At present there are two programs in use to compute the two-jet (plus beam 
remnant) cross sectiim in NLO using a JADEcj type of jet definition: PRO- 
JET and DIS JET EI There is now a new program, MEPJETp which can be 
used to compute jet cross sections using any infrared-safe definition. In addi- 
tion there is the promise of a generaUeurpose program for the NLO calculation 
of any infrared-safe quantity in DISu3 



Table 1: Predicted two-jet inclusive cross sections in DIS at LO and NLO, for four jet 
definition schemes and (at NLO) three different recombination schemes. 



Scheme 


LO 


NLO (E) 


NLO (EO) 


NLO (P) 


W 


1020 pb 


2082 pb 


1438 pb 


1315 pb 


JADE 


1020 pb 


1507 pb 


1387 pb 


1265 pb 


kx 


1067 pb 


1038 pb 


1014 pb 


944 pb 


cone 


1107 pb 


1203 pb 


1232 pb 


1208 pb 
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Results Ell on the inclusive two-jet cross section using MEPJET with a 
variety of jet definitions, for a particular set of acceptance cuts at HERA, are 
summarized in Table |l]. The jet-defining scheme that has been most widely 
used for NLO calculations is the so-called Vl^-schemepa in which two particles 
or jets are clustered if their invariant mass-squared Sy is less than j/cutW^^, 
where ?/cut is the jet resolution (j/cut = 0.02 here) and W is the hadronic cm. 
energy. Experimentally, the JADE clustering variable Mfj = 2EiEj{l — cos6i^ 
has been used in place of s^. A suggested improvement is the fc_L-schemepl 
in which fc^ = 2 inm{Ef , E^} {1 — cos 9ij), evaluated in the Breit frame of 
reference, is compared with a resolution scale (40 GeV^ here). Alternatively, 
one may use instead a cone type of jet definition, in which clustering occurs if 
(jii — TjjY -\- {(pi — 4>jY < 1, where 77 and (j) the pseudorapidity and azimuth 
in the HERA lab frame. In all these schemes the way of defining the combined 
momentum after clustering caxLalso be varied, leading to the so-called E, EO 
and P recombination schemesEj 

It can be seen from Table Q that in the W and JADE schemes there are 
large values of the if -factor {K = [NLO-LOJ/LO) and strong recombination 
scheme dependences, which are reduced in the kj_ and cone schemes. One 
worrying point is that the MEPJET results da not agree with those of PRO JET 
and DISJET for the W or JADE schemesEil This makes it very desirable to 
have another independent—calculation of these quantities, which is said to be 
coming in the near futuretJ 



5.2 Forward Jet Production 

The rate of production of fast (xjct S> x) forward jetsjii DIS at small x has 
been proposed as a good probe of BFKL dynamicsE^a The idea (Fig. ^ is 
that for jet transverse momenta Pt jet ^ there is little scope for ordinary 
DGLAP evolution, while the condition a;jct ^ x means that BFKL evolution 
in log X is enhanced. 

In order to show that there is a BFKL enhancement in forward jet pro- 
duction at small x, one should establish that the usual NLO calculation with 
DGLAP evolution is not able to explain the data. This is difficult because 
NLO processes like boson^luon fusion with an extra jet are dominant in the 
selected kinematic regionEd Thus NLO corrections are large and not well under 
control. 

The data of the HI collaborationEitogether with the results of a recent 



calculation based on BFKL dynamicsp3 are shown in Fig. 10 as a function 
of Bjorken x. The HI selection criteria for forward jets are a;jct > 0.025, 
0.5 < PT,jct/Q^ < 4, PTjct > 5 GeV and 6° < Ojct < 20°, where 9 jet is 
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Figure 9: Forward jet production as a probe of BFKL dynamics. 



the angle between the forward jet (defined by the cone algorithm) and the 
proton beam. The agreement between the BFKL calculation and the data is 
encouraging. The approximate matrix element calculation shown, performed 
by omitting the BFKL ladder in Fig. ^, falls well below the data. It will be 
important to check that this remains true when the full NLO matrix elements 
are included. For this purpose the new NLO programs discussed above, which 
allow the cone jet definition to be used, will be most valuable. 
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